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Introduction 


During  February  - March,  1975,  detailed  photometric  data  was  obtained 
along  +5°  of  the  geomagnetic  meridian  centered  at  the  base  of  the  ATS-5 
and  ATS-6  field  lines.  DMSP  satellite  images  were  also  available  for 
the  period. 

This  project  was  concerned  with: 

1.  Analysis  of  meridian-scanning  photometer  data  and  presentation  as  gray- 
scale "keogram"  plots. 

2.  Reduction  of  ATS-5  and  ATS-6  particle  data,  and  presentation  as  gray- 
scale electron  and  proton  "spectrograms." 

3.  Comparison  of  co-ordinated  photometer  and  DMSP  data  to  absolutely 
calibrate  the  DMSP  imaging  system. 

4.  Collection  of  all  relevant  ground-based  magnetograms. 

5.  Analysis  of  data  to  search  for  relationships  between  satellite  particle 
fluxes  and  energy,  and  auroral  intensity  at  the  base  of  the  same  field 
line.  In  particular,  to  examine  loss-cone  fluxes  as  measured  by  ATS-6 
for  this  purpose. 

6.  Publication  of  results. 
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Data  Analysis 


Table  1 gives  the  dates  for  which  data  were  analyzed.  These  were 
clear,  moonless  nights  for  which  both  ground-based  and  satellite  data 


were  complete. 


Table  1 


Day 

Date 

Keograms 

Spectrograms 

DMSP 

passes 

(1975) 

ATS-5 

ATS-6 

No. 

Time (IT) 

39 

Feb . 

8 

/ 

V 

/ 

/ 

4665 

0557 

4666 

0739 

40 

Feb. 

9 

/ 

/ 

/ 

4680 

0720 

43 

Feb. 

12 

/ 

/ 

/ 

4722 

0626 

4723 

0807 

44 

Feb. 

13 

/ 

/ 

/ 

4736 

0607 

4737 

0747 

45 

Feb. 

14 

/ 

X 

/ 

4751 

0729 

49 

Feb. 

18 

j 

X 

/ 

4807 

0614 

4808 

0756 

50 

Feb. 

19 

/ 

X 

/ 

4822 

0737 

64 

Mar. 

5 

/ 

X 

/ 

V 

5020 

0642 

5021 

0820 

66 

Mar . 

7 

/ 

X 

/ 

5049 

0744 

67 

Mar . 

8 

/ 

X 

/ 

5063 

0724 

68 

Mar. 

9 

j 

X 

/ 

5078 

0848 

71 

Mar . 

12 

/ 

X 

J 

5120 

0753 

72 

Mar. 

13 

/ 

X 

/ 

5134 

0734 

73 

Mar . 

14 

/ 

X 

/ 

5148 

0715 

74 

Mar. 

15 

j 

X 

/ 

5162 

0659 

1.  Photometric  Data 
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Analysis  of  photometric  data  was  a joint  computer  project  between 
Boston  College  and  Lockheed  Research  Labs  (under  subcontract).  These 
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these  it  was  decided  that  line  plots  of  various  satellite  measurements 
(energy  flux,  mean  energy)  were  necessary.  This  work  is  presently  in 
progress . 

It  is  hoped  these  additional  investigations  will  lead  to  new  pub- 
lishable results,  and  that  the  work  will  be  completed  over  the  next 
six  months.  As  this  contract  has  expired,  this  work  will  be  completed 
at  no  additional  cost.  However,  it  is  requested  that  publication  funds 
that  remain  unspent  in  the  contract  be  available  for  publication  charges 
if  an  additional  publication  results. 


* 


In  addition  to  the  calibration,  the  nature  of  particles  responsible 
for  the  "diffuse  aurora"  seen  on  the  low-latitude  edge  of  auroral  preci- 


pitation was  investigated.  It  was  found  at  times  the  observed  intensities 
could  be  explained  as  being  due  entirely  to  precipitating  protons.  At 
other  times,  electrons  contributed  up  to  2/3  of  the  total  excitation.  We 
believe  this  diffuse  aurora  represents  plasma  sheet  "drizzle,"  modified 
at  times  by  particle  drifts  after  substorm  rejection  events  (see  Appendix  I). 

4.  Magneto grams 

To  identify  substorm-associated  injection  events,  it  is  necessary  to 
collect  magnetograms  from  a number  of  stations  well-distributed  in  longi- 
tude. The  stations  used,  and  U.T.  of  local  midnight  for  each,  are  shown 
in  Fig.  5.  From  these  magnetograms,  substorms  were  identified  to  aid  in 
interpretation  of  the  relations  between  keograms  and  spectrograms. 

5.  Analysis  and  Publication 

Comparisons  of  these  data  have,  in  general,  confirmed  results  of  a 
similar  previous  study  (R.H.  Eather,  S.B.  Mende,  and  R.J.R.  Judge: 

Plasma  injection  at  synchronous  orbit  and  spatial  and  temporal  auroral 
morphology,  J.  Geophys . Res.,  81 , 2801,  1976).  Expected  new  insights 
with  this  new  data  set  have  not  materialized.  Much  of  the  problem  here 
is  that  ATS-6  did  not  sample  in  the  loss  cone  as  was  expected  when  this 
experiment  was  planned. 

An  experimenters  meeting  was  recently  held  at  UCSD  (late  June,  1978) 
with  Drs.  De  Forest,  Mcllwain,  Mende,  and  Eather.  It  was  agreed  that  it 
was  of  little  use  publishing  data  that  simply  confirmed  results  of  the 
previous  study  referenced  above.  However,  some  events  were  identified 
that  did  not  seem  to  fit  the  previous  picture.  To  further  investigate 
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tapes  had  numerous  errors  and  the  analysis  was  more  involved  and  took 


much  more  effort  and  time  than  expected. 

Keograms  were  prepared  at  three  wavelengths: 

a.  4278  N+  - a measure  of  total  energy  influx. 

2 

b.  6300  01  - a measure  of  low-energy  electron  influx. 

c.  4861  Hg  - a measure  of  proton  influx. 

Examples  of  a quiet  day  (March  7)  and  medium-activity  day  (Feb.  14) 
are  shown  as  Figs.  1 and  2. 

2.  Particle  Data 

Production  of  particle  spectrograms  was  done  under  subcontract  at 
University  of  California,  San  Diego.  Only  partial  data  was  available 
from  ATS-5.  The  greatest  disappointment  was  that  the  ATS-6  detectors, 
which  were  supposed  to  sample  the  particle  loss  cone,  were  always  at 
higher  pitch  angles.  Consequently  one  of  the  major  thrusts  of  this 
project,  that  of  comparing  average  energies  and  energy  fluxes  within  the 
loss  cone  (at  the  equator)  with  auroral  spectral  intensities,  could  not 
be  accomplished. 

Examples  of  the  ATS-6  spectrograms  for  the  "quiet"  and  "medium- 
activity"  days  (shown  in  Figs.  1 and  2)  are  shown  in  Figs.  3 and  4. 

3.  DMSP  Calibration 

A complete  and  detailed  calibration  of  the  DMSP  imaging  system  has 
been  achieved.  This  is  the  first  time  an  absolute  calibration  has  been 
obtained,  and  this  will  greatly  expand  the  usefulness  of  these  images  in 
auroral  and  magnetospheric  research.  A paper  reporting  these  results  has 
been  submitted  to  J.  Geophys.  Res,  and  is  attached  to  this  report  as 
Appendix  I. 
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DMSP  Calibration  and  the  Diffuse  Aurora 
or 

What  You  Always  Wanted  to  Know  About  DMSP 
But  Were  Afraid  to  Ask 


R.  H.  Eather 

Physics  Department 
Boston  College 
Chestnut  Hill,  MA  021A7 


Abstract 

Although  DMSP  satellite  data  is  widely  used,  there  has  been  no  reliable 
absolute  calibration.  Co-ordinated  data  with  ground-based  photometers 
allow  a calibration  curve  of  film  density  versus  4728  intensity  to 

be  derived.  The  DMSP  satellites  (5C  Series)  record  airglow,  and  can 
detect  auroral  forms  of  intensities  50  R of  4278  N9+-  It  is  estimated 
that  the  5D  Series  satellites  are  capable  of  detecting  auroras  with  25P 
of  4278  ,\'0+.  Comparison  of  "diffuse”  auroral  images  with  photometric 
data  show  that  this  form  of  aurora  is  excited  by  a range  of  proton  and 
electron  contributions,  varying  from  almost  pure  proton  excitation  to 
dominant  electron  excitation. 
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Introduction 

Since  1971,  satellite  imaging  of  the  aurora  has  been  used  to  monitor 
auroral  activity  on  a worldwide  basis.  The  Isis  2 spacecraft  has  provided 
spectral  imaging,  while  the  USAF  Air  Weather  Service  DAPP  and  DMSP  satellites 
have  provided  wide-band  imaging  (extending  from  the  visible  through  the  near 
infrared) . The  latter  images  are  currently  archived  and  are  available  from 
World  Data  Service  as  35  mm  microfilm  copies.  This  ready  access  has  resulted 
in  the  DMSP  photographs  becoming  an  important  resource  in  auroral  and  magneto- 
spheric  research. 

A major  problem  in  using  the  DMSP  images  has  been  the  lack  of  any  reliable 
calibration  data,  and  uncertainty  as  to  what  the  images  represent  in  terms  of 
the  traditionally  observed  auroral  emissions  in  the  visible  region.  Conse- 
quently the  DMSP  pictures  have  only  been  used  in  a qualitative  way,  which  is 
an  unfortunate  underutilization  of  a high-quality  system  that  has  the  potential 
for  accurate  quantitative  use.  The  purpose  of  this  paper  is  to  present  a 
detailed  calibration  for  one  of  a DMSP  satellite  (by  co-ordinated  comparisons 
with  ground-based  photometric  measurements) , and  illustrate  the  usefulness  of 
the  images  for  quantitative  investigations. 

The  DMSP  System 

The  DMSP  satellite  system  has  been  described  elsewhere  (Pike,  1975)  so  onlv 
relevant  operational  information  is  summarized  here.  The  system  consists  of 
two  satellites  in  sun-synchronous  polar  orbit  (altitude  830  km,  inclination 
98.75°,  period  101.5  mins.),  one  near  the  dawn-dusk  meridian  and  the  other  near 
the  noon-midnight  meridian.  Images  are  produced  by  a line-scanning  radiometer 
which  builds  up  pictures  by  repetitive  scanning  across  the  earth  along  the  orbital 
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track.  Resolution  is  3.7  1cm  (center  of  picture)  - 22  km  (edge  of 
picture)  for  the  earlier  Block  5B/5C  satellites  and  2.7  km  (center)  - 
5 km  (edge)  for  the  Block  5D  satellites  (operational  after  March  1977, 
see  Table  1).  "Resolution"  referred  to  here  is  not  simply  determined 
by  adjacent  scanline  spacing  (as  pixel  field  of  view  overlaps  slightly  at 
nadir,  and  more  at  the  edge  of  scan).  The  better  5D  edge  resolution  is 
accomplished  by  a step-function  decrease  of  field  of  view  and  variation 
of  the  angular  sampling  rate  to  compensate  for  increased  ground  area 
in  the  field  of  view  with  increasing  slant  range. 

The  following  discussion  refers  to  the  Block  5B/5C  satellites, 

A number  of  system  changes  were  incorporated  in  the  later  Block  5D 
satellites,  and  these  are  discussed  in  a later  section. 

The  "visual  channel"  used  for  auroral  images  was  sensitive  over 
a band  from  4,000  - 11,000  8 for  the  53/5C  satellites,  with  peak 
response  at  a wavelength  of  8000  8 (see  Figure  1).  Each  detect*  r had 
three  possible  gain  states  (selection  depending  on  overall  light  1 ev <. 1 ) , 
and  for  auroral  images  the  highest  gain  was  used  (except  in  strong 
moonlight  when  the  medium  gain  state  may  have  been  used).  The  earlier 
satellites  (Block  5B)  only  had  linear  amplifiers,  whereas  some  later 
satellites  (Block  5C)  have  both  linear  and  logarithmic  amplifiers  (see 
Figure  2).  Fiducial  markings  on  the  photographs  show  the  gain  state 
of  the  amplifier,  but  do  not  indicate  whether  the  amplifier  used  is  in 
linear  or  logarithmic  mode.  This  information  is  apparently  no  longer 
available  (even  from  the  Air  Weather  Service),  but  apart  from  some 


initial  testing  procedures,  it  appears  all  the  Block  5C  satellites 


were/are  in  the  logarithmic  mode  (and  highest  gain  state)  for  nighttime- 
auroral  imaging.  Caution  should  be  exercised,  however,  as  there  have 
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been  periods  in  recent  years  when  the  detector  was  in  the  linear  mode 
(because  of  system  testing  and  because  of  operational  errors). 

For  Block  5B/5C  satellite,  d.c.  voltage  levels  from  the  detector 
were  telemetered  to  the  ground,  where  they  were  digitalized  into  64 
levels,  and  stored  on  magnetic  tape.  (For  the  Block  5D  satellites,  128 
digital  levels  are  telemetered  directly  from  the  satellite,  but  the 
least  significant  bit  is  dropped  before  plotting.)  These  64  levels 
define  the  gray  scale  levels  for  the  resulting  transparencies.  The 
digital  data  on  tape  is  only  kept  for  a few  days  to  a week  before  tapes 
are  recycled. 

The  grav-scale  picture  plots  are  made  on  any  of  a number  of  plotters 
(three  are  currently  used).  Beginning  April,  1975  (just  after  the  data 
reported  here  were  obtained),  step-wedge  calibration  plots  were  provided 
daily  for  each  plotter  used  on  that  day.  The  plotter  used  for  each  DMSP 
picture  is  indicated  by  the  number  of  small  triangles  that  appear  (after 
June,  1976)  between  timing  information  along  tiie  edge  of  the  image.  To 
give  an  indication  of  the  variability  between  plotters,  gray-scale 
calibrations  for  three  plotters  used  on  April  18/19,  1978  were  used  to 
obtain  microdensitometer  plots  shown  in  Figure  3.  Also  shown  is  the 
nominal  gray  scale  calibration;  plotters  are  supposed  to  be  adjusted  to 
this  curve,  with  an  allowed  tolerance  of  +0.05D.  However,  it  has  been 


found  that  a density  curve  which  progresses  towards  white  faster  than 
the  nominal  linear  curve  produces  a much  "prettier"  picture  (clouds 
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appear  whiter  rather  than  gray)  so  the  plotters  are  often  aligned  to 
produce  this  characteristic.  Figure  3 indicates  the  CDM4  plotter  on  this 
day  was  adjusted  to  the  nominal  linear  curve,  whereas  the  CDM2  and  CDM3 
plotters  were  adjusted  for  "viewer  preference". 

It  was  also  found  that  the  gray  scale  calibration  curves  vary  some- 
what across  the  image.  This  variation  results  from  the  state  of  adjust- 
ment of  plotter  brightness  control  in  the  rectification  (linearization 
of  geocentric  angle  along  image  scanline)  alogrithim,  from  focus  varia- 
tions, changing  optics  path,  and  non-uniform  processing.  Consequently, 
for  the  most  accurate  work,  calibration  curves  should  probably  be 
derived  from  the  same  relative  location  as  the  data  region  of  interest. 

When  the  original  8"  positives  are  made,  they  are  copied  onto  35  mm 
positive  film  for  the  World  Data  Center,  and  the  originals  then  retained 
by  Air  Force  Geophysics  Labs.  Up  until  June,  1978,  WDC  treated  this 
35  mm  copy  as  their  master,  from  which  they  made  a positive  copy.  This 
positive  copy  was  then  used  to  make  copies  when  requested  by  users. 

Thus  users  obtaining  35  mm  microfilm  from  WDC  received  4th  generation 
film,  with  consequent  increase  in  contrast  and  loss  of  resolution.  At 
times  weak  auroras  clearly  visible  on  original  8"  positives  were 
completely  lost  (see  for  example  Figure  4a).  To  improve  quality,  as  of 
June,  1978,  WDC  makes  user  copies  from  their  35  mm  master,  so  the  user 
now  receives  3rd  generation  film. 

The  degree  of  contrast  enhancement  because  of  multiple  copying  is 
difficult  to  quantify,  as  various  film  types  and  processing  laboratories 
have  been  used  over  the  years.  Figure  4a  illustrates  the  problem  is 
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serious  at  times.  However,  a comparison  done  in  this  study  of  micro- 
densitometer scans  of  an  original  8"  positive  and  a 4th  generation  35  mm 
positive  obtained  from  V.'DC  (Figure  4b)  shows  surprisingly  good  reproduction. 
The  35  mm  copy  is  clearly  of  little  use  for  quantitative  work,  but  most 
auroral  features  have  been  preserved. 

Previous  Calibrations 

One  of  the  earliest  reports  of  the  DMSP  pictures  (Rogers  et.  a]., 

1974)  stated  that  the  minimum  detectable  signal  was  8.5  x 30  10  watt 

-2-1  -Q  _o_i 

cm  ster  , and  that  the  system  saturated  at  8.5  x 10  watt  cm  “ ster 
The  lower  number  has  no  analytical  basis  and  rests  on  the  subjective 
assumption  that  the  lowest  107  of  the  scale  at  high  gain  would  be  "very 
noisy"  and  so  of  "little  value".  These  spectral  integrated  fluxes  are 
for  a solar  spectrum  (0.4  - 1.1„)  and  correspond  to  about  43  kR  and 
430kR  respectively  (along  the  line  of  sight  of  the  instrument).  These 
numbers  refer  to  the  #5528  Block  5B  satellite,  which  only  had  linear 
amplifiers,  and  so  covered  only  one  order  of  magnitude  dynamic  range 
(Figure  2).  Although  "typical"  for  the  5B/5C  satellites,  there  was  sotth 
variation  between  the  various  satellites. 

Berkey  and  Kamide  (1976)  used  predicted  intensities  of  aurora] 
emissions  (from  Vallance  Jones,  1971)  within  the  pass-band  of  the  DMSP 
detector  to  determine  that  the  threshold  sensitivity  corresponded  to 
2-3  kR  of  the  5577  01  emission  and  saturation  corresponded  to  30  kR. 

A number  of  nearly  coincident  passes  of  the  DMSP  and  Isis  2 satellites 
allowed  them  to  attempt  intercalibration,  from  which  they  concluded  the 
DMSP  radiometer  threshold  corresponded  to  1.7  - 3.3  kR  of  3914  N0+ 
emission.  They  warned  that  there  is  considerable  emission  below  this 
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threshold  and  that  investigators  using  DMSP  imagery  should  be  cognizant 
of  that  fact. 

These  numbers  suggest  that  the  DMSP  detector  is  not  as  sensitive 
as  all-sky  cameras,  which  typically  give  useable  images  corresponding 
to  '.75  kR  of  5577  01.  In  the  mid-1970's,  it  became  apparent  that  DMSP 
images  were,  if  anything,  somewhat  more  sensitive  than  all-sky  cameras. 

Such  conclusions  came  from  qualitative  comparison  of  photographs,  and 
were  transmitted  by  word  of  mouth  rather  than  by  formal  publication. 

This  increased  sensitivity  of  DMSP  clearl  resulted  from  the  new 

logarithmic  amplifiers.  Reference  to  Figure  2 shows  that  the  8.5  x 10 

-9  -2  -1 

• 8.5  x 10  watts  cm  ster  operational  range  of  the  linear  detector 

results  in  a channel  output  voltage  of  0.6  *7.2  volts.  This  same 

-10  -9 

output  voltage  range  corresponds  to  1.6  x 10  -*  8.5  x 10  watts 

-2  -1 

cm  ster  for  the  logarithmic  amplifier,  representing  an  increase  in 
sensitivity  of  a factor  of  5,  an  increase  in  dynamic  range  to~50,  and 
the  same  saturation  level.  Thus  according  to  Berkev  and  Kamides  (.1976) 
numbers,  we  would  expect  a threshold  level  of  300-600  R at  3914  N.,+  , 
and  saturation  at  25  kR  (uncorrected  for  albedo). 

Mende  and  Eather  (1976)  reported  an  intercalibration  of  DMSP  (“8531 
with  logarithmic  amplifier)  with  ground  based  photometers.  They  con- 
cluded that  DMSP  gave  a useable  image  when  auroral  intensity  reaches 
150  R at  4278  N^+  (corresponds  to  450  R at  3914  N?+) , but  that  the 
picture  saturated  at  '500  R 4278  N9+.  It  will  be  shown  in  this  paper 
that  DMSP  sensitivity  is  in  fact  as  low  as  50  R at  4278  N^+,  and  saturation 
occurs  at  '2.5  kR.  The  discrepancy  witli  the  previous  work  results  from 
the  present  use  of  original  8"  formal  DMSP  positives,  rather  than  the  4th 
generation  copy  prints  obtained  from  World  Data  Center  microfilm  for  the 
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earlier  study. 

Airglow  and  Continuum 

With  the  wide  spectral  band  detectors  used  on  DMSP,  the  integrated 
airglow  and  continuum  emissions  is  appreciable.  The  principal  airglow 
emissions  are  the  atomic  oxygen  lines,  the  0^  atmospheric  system,  and 
the  hydroxyl  system.  Expected  intensities  (Roach  and  Gordon,  1973)  are 
shown  in  Figure  lb,  together  with  the  detector  response  curve.  Integra- 
tion over  the  response  curve  gives  about  18  kR.  To  this  should  be  added 
that  contribution  from  reflected  starlight,  zodiacal  light  and  gegen- 
schein,  which  might  be  expected  to  add  another  2-3  kR  (Roach  and  Gordon, 
1973).  Assuming  an  albedo  of  1.0  for  the  snow  covered  arctic,  and 
allowing  for  atmospheric  transmission  effects  (Xawrocki  and  Papa,  1963), 
we  estimate  that  DMSP  will  see  about  30-35  kR.  This  total  "background" 
emission  will  trigger  the  first  307  of  the  dynamic  ranee  of  the  system 
(see  Figure  2)  and  so  result  in  an  overall  background  density  on  the 
film  of  '1.15  D (See  Figure  3). 

The  above  comments  apply  on  moonless  nights.  Background  levels 
will,  of  course,  rise  dramatically  under  lunar  illumination  and  result 
in  considerable  loss  in  ability  to  detect  weak  auroral  emissions. 

At  this  point  it  is  of  interest  to  calculate  the  auroral  intensitv 
needed  to  equal  this  airglow/continuum  contribution.  Reference  to 
Figure  la  gives  that  a total  aurora]  intensity  of  33  kR  (integrated 
over  detector  response  and  with  some  albedo  and  transmission  assumptions 
as  above)  corresponds  to  a 5577  01  intensity  of  670  R.  Reference  to 


Figures  2 and  3 shows  that  if  this  is  added  to  the  airglow/continuum 
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the  DMSP  image  density  should  change  by  about  0.2D  to  0.95D.  Weak 
auroral  images  will  appear  as  low-contrast  enhancements  (diffuse 
borders)  on  the  background  density,  which  is  typically  quite  grainy. 

A human  observer  should  be  able  to  readily  identify  a 0.05D  change 
under  such  conditions,  which  corresponds  to  170  R 5577  01 . This 
should  represent  an  estimate  of  the  ultimate  sensitivity  of  the  DMSP 
system  i.e.,  an  aurora  of  '170  R 5577  01  should  be  just  detectable. 
(This  corresponds  to  100  R 3914  N?+  or  35  R 4278  N +.) 

A four-channel  meridian-scanning  til t ing-f i 1 t er  photometer  system 
was  operated  from  Gillam,  Manitoba  (6b.2°N,  323. 4°E  geomagnetic)  in 
February-March , 1975  for  co-ordinated  measurements  under  the  ATS-5  and 
ATS-6  satellites.  During  this  period,  there  were  20  cases  of  simul- 
taneous DMSP  and  photometer  data  on  clear,  moonless  nights,  with  Gillam 
appearing  near  the  central  part  of  the  DMSP  picture  and  auroras  in  the 
field  of  view. 

The  photometer  scanned  th (V  from  the  zenith,  with  a nc th-south 
scan  taking  90  secs.  The  time  of  DMSP  passage  across  the  Gillam 
meridian  was  determined,  and  the  photometer  scans  immediately  before, 
during,  and  after  this  time  were  used  to  compare  with  m : crodensitomet or 
plots  from  the  original  DMSP  8"  positives.  The  mic rodensi tometer  scan 
was  made  along  the  projection  of  the  photometer  meridian  scan  at  the 
100  km  level  on  the  DMSP  picutre.  (Gridding  procedures  are  explained 
by  Pike,  1975.)  Microdensitometer  spot  size  was  adjusted  to  roughlv 
match  the  3°  field  of  view  of  the  photometers.  (In  most  cases  it  was 
possible  to  identify  the  town  lights  of  Gillam  on  the  DMSP  picture. 


and  the  image  of  these  lights  gave  a corresponding  peak  on  the  micro- 


11 


densitometer  scan.)  Photometric  ratios  were  used  to  derive  heights  of 
emission,  and  hence  latitude,  as  in  our  standard  keogram  program  (Eather 
et.  al . , 1975).  Calculated  field-aligned  photometric  intensity  can  then 
be  plotted  as  a function  of  latitude. 

A logarithmic  density  range  of  0.05  -»  2. ID  was  necessary  to  span 
the  transmission  range  from  the  clearest  part  of  the  DMSP  positive  (fog 
level)  to  fully  saturated  black.  The  saturated  black  was  down  the  sides 
of  the  film  (various  fiducial  marks  etc.),  an  the  "blackest"  part  of 
the  DMSP  images  corresponded  to  a density  of  1.2  I):  this  density  level 
corresponds  to  the  integrated  airglow  signal. 

"Weak"  Auroral  Arcs:  An  example  of  co-ordinated  data  for  very  weak 
auroral  conditions  is  shown  in  Figure  5.  The  region  of  the  photometer 
scan  is  marked  on  the  DMSP  picture,  which  is  a contact  print  from  the 
original  positive.  The  diffuse  aurora  in  the  bottom  part  of  the  picture 
is  stable  and  does  not  change  during  the  three  photometer  scans.  The 
structured  arc  at  the  top  of  the  picture  is  clearly  moving  and  transient, 
as  it  only  shows  in  the  photometer  scan  before  the  co-ordination  time, 
and  then  was  at  a lower  latitude  than  when  detected  hv  DMSP. 

The  overall  background  density  of  the  DMSP  positive  is  about  1 . 2D 
(and  this  was  true  for  all  positives  examined).  The  diffuse  arc  image 
appears  above  this,  and  it  corresponds  well  with  the  photometer  data. 

4278  N’2+  intensities  exceeding  40  R are  seen  to  give  a DMSP  densitv 
above  the  1 . 2D  background.  Peak  4278  N0+  intensity  of  1 20  R corresponds 
to  a film  density  of  about  .950. 

Also  plotted  in  Figure  5 is  the  H intensity.  The  ratio  of 
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4278  N2+/H  is  '1.5,  which  implies  that  all  of  the  observed  diffuse  aurora 
can  be  excited  by  precipitating  protons  (Eather,  1968).  The  narrow 
structured  arc  on  the  northern  edge  of  the  diffuse  aurora  is  too  narrow 
to  be  excited  by  protons  (because  of  charge-exchange  spreading),  so  that 
arc,  together  with  the  transient  arc  to  the  north,  must  be  primarily 
electron  excited. 

Three  more  examples  of  weak  arcs  are  given  in  Figure  6a,  b,  c.  In 
Figure  6a,  the  double  structure  of  the  diffuse  arc  is  not  resolved  by 
the  photometer  keogram  program  (which  is  understandable  if  one  considers 
the  vertical  distribution  of  emissions  and  geometrv  of  the  arcs  with 
respect  to  the  station).  The  H plot  and  4278  N0+/H  ratio  of  4.  show, 
that  both  electrons  and  protons  contribute  to  diffuse  aurora  excitation, 
though  the  region  of  hydrogen  emission  is  broader.  In  Figure  6b,  the 
northernmost  intense  structure  is  clearly  transient;  the  patchy  aurora 
over  Oillam  is  also  undergoing  short  term  changes;  the  diffuse  aurora 
to  the  south  is  stable,  and  the  ratio  of  4278  N'7+/H  ( 4)  indicates 

protons  and  electrons  contribute  about  equally  to  its  excitation,  with 
the  hydrogen  emissions  broader  in  latitude.  In  Figure  6c.  the  strong 
northern  arc  is  transient  and  must  have  only  existed  by  -90  seconds, 
just  as  the  satellite  was  passing  over  that  location.  The  stable 
diffuse  aurora  seems  to  be  largely  excited  by  electrons  on  its  northern 
half  and  protons  on  the  southern  half. 

"Medium"  Auroral  Arcs:  Figure  7a,  b,  c gives  examples  of  more 
intense  auroral  forms.  In  Figure  7a,  there  is  an  overall  background 
level  of  180  R of  4278  N0+  which  gives  a film  density  of  -0.83P. 
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About  1/3  of  this  seems  to  be  proton  excited,  and  the  remainder  must  be 
widespread  low  intensity  electron  precipitation.  Superimposed  on  this 
to  the  south  of  Gillam  is  a double-structured  diffuse  form  excited  by  a 
combination  of  electrons  and  protons,  and  temporal  changes  are  evident 
in  theelectron  component.  In  the  example  of  Figure  7b,  most  of  the 
-300  R 4278  N?  across  the  picture  can  be  attributed  to  the  broad  diffuse 
proton  precipitation,  giving  H -100  R.  The  electron-excited  structures 
that  superimpose  vary  in  time,  so  detailed  correlations  are  difficult  to 
make,  though  the  arc  almost  overhead  at  Gillam  is  identifiable  in  the 
DMSP  and  photometer  trace  corresponding  to  the  satellite  being  overhead. 
4278  intensity  of  1000R  gives  a film  density  of  .370.  Figure  7c 
shows  a localized  temporal  structure  in  the  north  (not  seen  by  the  photo- 
meter), a changing  but  well  correlated  structure  in  the  center  of  the 
picture,  and  a stable  diffuse  aurora  to  the  south  which  seems  dominantly 
proton  excited  but  with  some  superimposed  electron  structure. 

"Strong"  Auroral  Arcs:  Figure  8a,  b,  c shows  examples  of  auroras 
that  are  intense  enough  to  saturate  the  DMSP  film.  In  Figure  10a,  the 
DMSP  is  just  on  the  verge  of  saturation,  and  this  is  seen  to  correspond 
to  a 4278  N?+  intensity  of  2.75  kR.  The  strong  diffuse  aurora  in  the 
south  ( 1 30  R Hp,  800  R 4278  N^"1")  is  excited  by  both  protons  and  electrons. 
Figures  8b  and  8 c show  similar  examples,  with  some  DMSP  peaks  being 
transient  auroras  that  are  missed  by  the  photometer. 

Calibration  Curve 

The  20  cases  of  co-ordination  (of  which  10  have  been  presented  here) 
were  examined  and  all  cases  of  definite  correlated  structures  used  to  plot 
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the  calibration  curve  of  Figure  9.  As  previously  mentioned,  no  gray-scale 
calibration  data  was  available  at  this  time.  Fortunately  all  DMSP  pictures 
in  this  study  were  generated  on  the  same  plotter,  so  there  is  not  the 
uncertainty  of  different  calibrations  for  different  plotters  (Figure  3). 
However,  part  of  the  scatter  in  Figure  9 could  result  from  a changing 
calibration  of  the  plotter  with  time,  as  these  data  cover  a six  week 
period.  The  continuous-drive  wet  process  used  does  not  receive  good 
quality  control  and  maintenance,  and  chemical  residues  in  the  system  are 
common,  sometimes  leading  to  solid  matter  on  the  dried  films. 

It  may  be  seen  that  a consistent  calibration  relationship  is 
apparent.  We  conclude  that  absolute  intensities  may  be  derived  from  tin 
original  DMSP  positives,  within  an  accuracy  of  1501.  Minimum  detectable 


intensity  f- 


tellite  was  50  R of  4278  M,,  , with  saturation  of 


2.5  kR.  The  DMa.  ^ta  has  not  been  corrected  for  albedo,  which  is 
expected  to  be  1.0  for  the  uniform  snow  covered  surface  of  northern 
Manitoba.  This  effectively  increases  the  absolute  sensitivity  of  DMSP, 
but  as  albedo  affects  the  background  airglow  signal  also,  it  does  not 
change  the  signal/noise  ratio  that  determines  the  minimum  detectable 
auroral  intensity.  We  argued  above  (in  the  section  of  Previous  Calibra- 
tions) that  the  expected  threshold  for  the  logarithmic  detectors  should 
be  300-600  R at  3914  which  corresponds  to  100-200  R at  4278  N9+. 

Given  the  effect  of  albedo,  this  corresponds  to  50-100  R as  measured  for 
the  same  aurora  from  the  ground.  The  agreement  of  this  estimate  with  the 


experimentally  determined  50  R is  better  than  anyone  had  any  right  to 


expect . 
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Note  too  that  there  is  a van  Rhijn  effect  (the  increase  in  optical 
thickness  because  of  increased  obliquity  of  the  view  angle)  in  the  DMSP 
photographs.  Depending  on  the  height  distribution  of  the  aurora  and 
the  geometry  of  the  situation,  this  can  give  a factor  of  2 enhancement 
near  the  east  and  west  edges  of  the  picture.  This  was  not  allowed  for 
in  this  analysis  as  in  most  cases  Gillam  was  nearer  the  center  of  the 
picture.  (The  effect  is,  however,  allowed  for  in  the  program  that 
analyzes  the  photometer  data  - this  is  necessary  because  of  high 
obliquity  in  ground  based  measurements  for  aurora  still  near  the  central 
part  of  the  DMSP  pictures.) 

Warning  Concerning  Temporal  Fluctuations 

I 

] Ke  have  mentioned  a number  of  cases  above  where  auroral  forms  or, 

« » 

the  DMSP  pictures  are  not  recorded  on  the  photometer  scans,  implying 

\ 

that  they  existed  for  '90  seconds  (the  meridian  scan  time  of  the  photo- 

i 

meter),  so  appeared  briefly  at  their  recorded  position  when  DMSP  was 
. overhead  at  that  location.  In  most  cases  both  temporal  and  spatial 

effects  are  probably  important. 

Usually  examination  of  the  shape  of  the  form  indicates  it  was 
probably  a transient  brightening  (see  for  example  the  northern  forms  in 
Figures  Ab,  5,  6)  but  sometimes  what  looks  like  a stable  arc  form  on 
the  photograph  was  in  fact  a transient  phenomena.  Such  was  the  case 
published  by  Mende  and  Eather  (1976)  and  repeated  here  as  Figure  10. 

The  well  defined  arc  just  south  of  the  zenith  at  Gillam  in  the  picture 
actually  formed  north  of  Gillam,  quickly  moved  south,  then  faded,  all 


within  a 3 minute  time  interval.  There  were  no  other  discrete  auroras 
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recorded  by  the  photometer  within  t20  minutes  of  this  event.  This 
example  is  presented  simply  to  illustrate  the  possible  pitfalls  in 
assuming  longevity  and  stability  of  discrete  auroras  recorded  by  DMSP. 

Subsequent  Satellites 

The  calibration  curve  presented  here  is  for  the  5C/8531  satellite-. 

The  other  5C  satellites  (/•  9531  and  #10533)  also  used  radiometer  detectors, 
and  should  have  very  similar  calibration  curves.  The  5D  satellites  havi 
photomultiplier  detectors  with  less  red  response  than  the  radiometer 
(silicon  diode)  detectors.  Although  the  quantum  efficiency  of  photo- 
multipliers (peak  -'•25/  ) is  much  lower  than  for  silicon  diodes,  the 
signal/noise  characteristics  at  low  light  levels  are  far  superior. 

Normalized  response  curves  for  some  5D  detector  channels  (including 
optics),  are  shown  in  Figure  11a.  FI  to  F3  used  modified  extended  S-20 
photocathodes,  whereas  the  F4  and  F5  instruments  use  cesiated  gallium- 
arsenide  photocathodes  for  extended  red  response,  and  a yellow  filter 
to  limit  blue  response.  Peak  quantum  efficiencies  are  ~22  for  tri- 
alkali tubes,  and  — -1 52  for  GaAs  tubes. 

The  calibration  curve  for  the  F3  detector  shown  in  Figure  lib  illu- 
strates the  characteristics  of  the  5D  system.  The  photomultiplier 

detector  is  switched  in  line  when  radiance  falls  below  about  5 x 10  ^ 

-2  -1 

watts  cm  ster  , and  then  switched  to  a higher  gain  state  at  about 
-7  -2  -1 

4 x 10  watts  cm  ster  . (Radiance  is  a spectral  integration  over 
the  solar  spectrum  between  0.4  and  l.lq.)  Plotted  in  Figure  lib  is 


t 


r 


the  Variable  Digital  Gain  Amplifier  (VDGA)  gain,  which  lias  a 64db  dynamic 


range  (dB  = 20  x log  (ratio)).  For  a given  VDGA  gain,  the  plot  gives  the 

radiance  that  results  in  saturation  of  the  amplifier. 

Normally  for  auroral  images,  the  PMT  detector  is  in  a logarithmic 

mode,  with  a dynamic  range  of  40db.  Thus  at  face  value  Figure  lib  shows 

that  when  the  VDGA  amplifier  gain  is  near  its  maximum  of  fS4db,  the 

-9  -2  -1 

detectable  input  would  be  40db  down  on  ~5  x 10  watts  cm  ster 
-11  -2  -1 

i.e.  -5  x 10  watts  cm  ster  . However,  detector  noise  has  been 

measured  in  orbit  with  a known  intensity  LED  source,  which  gives  a S/N 

-10  2 

ratio  of  1.0  for  a radiance  of  ^5  x 10  watts/cirT  ster.  Expected 
airglow  levels  are  about  1.5  times  this  radiance,  and  so  effectively 
set  the  low  limit  for  auroral  detectability. 

To  estimate  the  relative  sensitivity  of  the  5D  and  5C  series  instru- 
ments to  an  auroral  spectrum,  we  folded  the  detector  spectral  responses 
into  the  predicted  auroral  and  airglow/continuum  spectra  shown  in  Figs, 
la,  lb,  and  Table  1 shows  the  results. 


Satellite  Number 


5C/8531  I 
5C/9532  / 
5C/1053  ) 

5D/12535  (FI)  | 
5D/13536  (F2)  / 
5D/14537  (F3)  ) 

5D/15538  ) , 


Table  1 


Spectral  Integrated  Flux  (kR) 


IBC  2 Aurora 


Airglow 


Contribution 
of  S"?  01 


5D/15PJS  ) , . . ) 

5D/16539  n0t 

' 1 aunr.ned  ) " 


— , launched)' 

Table  1 shows  that  the  reduction  in  auroral  response  of  the  5D  detector 
(because  of  its  decreased  spectral  range)  is  a factor  of  2 - 2.5.  However 
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the  relative  contribution  of  aurora  to  airglow  is  better  for  the  5D  detec- 
tor. If  we  define  the  "ultimate  auroral  sensitivity"  of  the  detector  as 
its  ability  to  detect  low  auroral  intensities  above  airglow,  we  estimate' 
the  current  5D  PMT  detectors  should  be  perhaps  twice  as  good  as  the  5C 
detector,  and  so  should  be  capable  of  detecting  aurora  characterized  by 
~25R  of  4278  • Future  5D  detectors  with  increased  red  response  should 

be  capable  of  detecting  aurora  characterized  by  ~35-40R  4278  N0+. 

Unfortunately  some  other  features  of  the  5D  instruments  make  it 
less  desirable  as  an  auroral  imager.  There  are  three  modes  of  gain  contr. 

(i)  Along  scan:  gain  is  adjusted  for  each  pixel  from  a program  that 

calculates  solar  elevation  wrt  viewing  region.  At 
quarter  moon  and  above,  expected  lun,.r  illumination 
as  a function  of  solar  elevation  is  calculated,  and 
used  to  control  gain. 

(ii)  Along  track:  gain  is  optimized  for  solar/lunar  elevation  each  scan 

at  subsatellite  point. 

(iii)  Preset:  gain  is  constant  for  picture. 

Modes  (i)  and  (ii)  are  to  optimize  cloud  pictures  and  will  normally  be 
operative  for  solar  elevations  above  -28°,  though  may  also  be  in  use 
for  some  auroral  images  (especially  if  there  is  lunar  illumination).  Such 
automatic  gain  changes  (in  increments  of  1/8  db)  would  make  it  virtually 
impossible  to  derive  quantitative  auroral  intensities.  But  normally  for 
new  moon  conditions,  the  detector  should  be  in  mode  (iii)  for  auroral 
images,  and  there  are  64  possible  preset  gains,  adjustable  from  0 to  63db 
in  ldb  steps.  The  gain  mode  and  gain  state  of  the  detector  is  coded  dovm 
the  side  of  the  DMSP  positives,  with  h;db  resolution  (see  Figure  12). 

Another  problem  arises  from  switching  of  the  pixel  field  of  view 
on  the  5D  satellites.  Because  of  increasing  distance  to  the  earth  at  the 
edges  of  the  pictures,  the  resolution  at  constant  field  of  view  would  be 


a 
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down  by  a factor  of  6 at  the  picture  edge  (as  were  the  5C  pictures).  The 
5D  instrument  changes  to  a smaller  field  of  view  about  1/3  of  the  way 
out  from  the  center  of  the  picture,  and  reduces  the  variation  in  pixel 
resolution  to  a factor  of  2.  When  field  of  view  is  decreased,  less  light 
enters  the  detector  so  amplification  is  simultaneously  adjusted  to  compen- 
sate. This  compensation  has  never  been  particularly  successful  and  as  a 
result,  the  center  part  of  auroral  pictures  is  effectively  at  a different 
gain  than  the  edges.  This  step  change  is  very  obvious  (see  Figure  13) 


an.:  would  have  to  be  allowed  for  in  quantitative  work. 
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Conclusions 

Tracking  down  the  chain  of  events  and  calibration  variables  between 
the  DMSP  detector  and  the  film  supplied  to  users  has  been  interesting 
detective  work.  There  is  no  central  documentation  source  and  much  of  the 
information  has  been  uncovered  by  word  of  mouth.  Perhaps  this  situation 
is  not  surprising;  the  DMSP  system  is  a real-time  operating  Air  Force 
system  which,  besides  supplying  meteorological  cloud  data,  is  used  to 
provide  auroral  positional  data  to  an  ionospheric  prediction  program  for 
real-time  communications  applications.  The  system  was  not  designed  for 
quantitative  auroral  reasearch,  and  the  Air  Force  has  been  most  co-opera- 
tive  in  making  auroral  data  available  to  the  scientific  community. 

The  so-called  "diffuse  aurora"  was  originally  identified  and  cate- 
gorized from  ISIS  pictures  (Lui  et.  al.,  1973).  From  the  characteristics  | 

and  location  of  this  diffuse  belt,  it  seemed  that  it  was  the  region  of 
proton  aurora,  and  insistence  on  a new  term  to  describe  this  region 
(simply  from  its  appearance  on  the  DMSP  photographs)  has  led  to  considerable 
confusion.  Nagata  et.  al . (1975)  demonstrated  that  the  diffuse  auroral 

belt  always  coincided  with  the  H emission  belt,  and  that  the  Correspon- 
ds 

dence  was  "perfect  and  without  exception."  The  present  study  confirms 
that  conclusion,  but  the  quantitative  data  shows  that  protons  are  not 
always  the  dominant  energy  source  in  this  region.  We  believe  this 
"diffuse  aurora"  is  a combination  of  steady  plasma-sheet  drizzle  of 


protons  and  electrons,  modified  after  injection  events  by  particle  loss 
from  the  reconstituted  plasma  sheet  of  freshly  injected,  azimuthally 
drifting  protons  and  electrons  (Eather  et.  al.,  1975).  Proton  drift 


paths  after  injections  can  lead  to  regions  of  essentially  pure  proton 
precipitation.  Consequently  the  proton/electron  ratio  of  precipitating 
particles  in  proton  auroras  (or  diffuse  auroras)  is  variable  and  depends 
on  the  recent  time  history  of  substorm  injections. 

The  DMSP  detector  is  a high  quality  instrument,  and  can  be  used 
profitably  in  a quantitative  fashion  if  users  are  aware  of  possible  pit- 
falls.  If  users  know  in  advance  of  requirements  for  data  at  special  times 
a request  might  be  made  to  assure  a desired  gain  mode  and  gain  setting, 
and  to  preserve  the  digital  tapes  for  superior  quantitative  data.  Compute 
processing  of  such  data,  or  of  digitized  data  regenerated  from  DMSP  posi- 
tives, could  be  used  to  determine  total  energy  inputs,  or  to  derive  inten- 
sity contour  maps  for  quantitative  comparison  with  other  measurements. 
Thus,  reliable  absolute  calibrations  can  greatly  expand  the  usefulness  of 
these  photographs  in  auroral  and  magnetosphc-r ic  research. 
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Fig.  la. 


Fig.  lb. 


Fig.  2. 


Fig.  3. 


Fig.  4a. 


Fig.  4b. 


Figure  Captions 


Intensities  of  the  main  auroral  emissions,  for  IBC2  Aurora  (5577  01  = 
10  kR) , from  Vallance  Jones  (1974).  Also  shown  is  the  spectral 
response  of  the  DMSP  radiometer  used  on  Block  5B/5C  satellites.  The 
histogram  shows  the  fraction  of  the  total  auroral  response  in  1000  X 
wavelength  intervals.  The  dashed  curve  is  a typical  response  charac- 
teristic (normalized)  for  all-sky  camera  (panchromatic)  film. 


Intensities  of  the  main  airglow  emissions,  from  Roach  and  Gordon 
(1973).  Also  shown  is  the  spectral  response  of  the  DMSP  radiometer 
used  on  Block  5B/5C.  satellites. 


Representative  response  curves  of  the  linear  and  logarithmic  ampli- 
fiers on  the  Block  5B/5C  DMSP  satellites.  The  radiance  is  spectrally 
integrated  over  a solar  spectrum. 


Film  density  versus  step  wedge  number  for  the  calibration  film  supplied 
with  DMSP  film  on  April  18/19,  1978.  The  CDM  number  defines  the 
plotter  used  to  generate  the  calibration  gray  scale. 


Contact  print  from  original  DMSP  8"  positive,  shows  weak  aurora  over 
Greenland  on  December  28,  1976.  On  the  right  is  an  enlargement  of 
the  copy  frame  provided  on  35  mm  microfilm  from  World  Data  Center. 


Comparison  of  microdensitometer  plots  from  original  DMSP  8"  positive. 


and  from  35  mm  microfilm  obtained  from  World  Data  Center,  for 
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February  13,  1975.  Also  shown  is  a contact  print  from  the  original, 
and  an  enlargement  from  the  microfilm. 


Fig.  5.  Contact  print  from  original  DMSP  8"  positive,  microdensitometer 

plot,  and  corresponding  ground  based  photometer  data,  for  Orbit 
4751  of  satellite  8531,  February  14,  1975.  The  meridian-scanning 
photometer  scans  along  the  direction  marked  on  the  DMSP  photograph, 
which  is  the  track  scanned  by  the  microdensitometer.  The  dashed 
peak  on  the  micro-densitometer  plot  is  the  town  lights  of  Gil  lam 
and  a nearby  dam  construction  site,  seen  as  two  dots  in  the  center 
of  the  meridian  scan  marked  on  the  DMSP  picture.  Three  photometer 

scans  at  4278  N + are  showm:  is  the  scan  just  before  DMSP  was 

overhead,  is  the  scan  during  which  DMSP  passed  overhead,  xxx 

is  the  scan  just  after  DMSP  was  overhead.  The  H region  is  shaded. 


Fig.  6a, 
b,c 


Further  examples  of  coordinated  data  for  "weak"  auroral  forms.  Sic 
legend  for  Figure  5. 


Fig.  7a, 
b , c 


Examples  of  coordinated  data  for  "medium"  intensity  auroral  f.  ms. 
See  legend  for  Figure  5. 


Fig.  8a, 
b,c 


Examples  of  coordinated  data  for  "strong"  auroral  forms.  Se> 
legend  for  Figure  5. 


Fig.  9.  Calibration  curve  derived  from  coordinated  data  where  auroral  forms 

are  clearly  correlated  on  the  DMSP  and  photometer  data. 


I 


Fig.  10  . 


DMSP  picture  for  February  13,  1975.  Also  shoum  arc  representative 
photometer  scans  before,  during,  and  after  the  satellite  overpass. 
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Fig.  11a. 


Fig.  lib. 


Fig.  12. 


; 

\ 

\ , . 

Fig.  13. 


I 


to  illustrate  the  temporal  nature  and  spatial  movement  of  the 
seemingly  steady  arc  seen  in  the  picture. 

Spectral  response  of  some  of  the  PMT  detectors  used  for  the  Block.  5!) 


satellites. 

Amplifier  gain  required  to  give  100"  output  (saturation)  for  the 
indicated  radiances,  for  the  F3  detector.  The  radiance  is  spectrally 
integrated  over  a solar  spectrum. 

Explanation  of  coding  of  time,  gain  mode,  gain  level,  and  gain 
trend  on  Block  5D  images.  This  coding  strip  appears  along  the 
right-hand  edge  of  the  images. 

Examples  of  discontinuities  introduced  by  switching  of  pixel  field 
of  view. 
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